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Abstract 

We investigate the equation of state, diffusion coefficient, and structural order of a family of 
spherically-symmetric potentials consisting of a hard core and a linear repulsive ramp. This generic 
potential has two characteristic length scales: the hard and soft core diameters. The family of 
potentials is generated by varying their ratio, A. We find negative thermal expansion (thermody- 
namic anomaly) and an increase of the diffusion coefficient upon isothermal compression (dynamic 
anomaly) for < A < 6/7. As in water, the regions where these anomalies occur are nested domes 
in the (T, p) or (T, P) planes, with the thermodynamic anomaly dome contained entirely within 
the dynamic anomaly dome. We calculate translational and orientational order parameters {t and 
Qq), and project equilibrium state points onto the {t,QQ) plane, or order map. The order map 
evolves from water- like behavior to hard-sphere- like behavior upon varying A between 4/7 and 6/7. 
Thus, we traverse the range of liquid behavior encompassed by hard spheres (A = 1) and water-like 
(A ~ 4/7) with a family of tunable spherically-symmetric potentials by simply varying the ratio of 
hard to soft-core diameters. Although dynamic and thermodynamic anomalies occur almost across 
the entire range < A < 1, water- like structural anomalies (i.e., decrease in both t and Qq upon 
compression and strictly correlated t and Qq in the anomalous region) occur only around A = 4/7. 
Water-like anomalies in structure, dynamics and thermodynamics arise solely due to the existence 
of two length scales, orientation-dependent interactions being absent by design. 
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I. INTRODUCTION 



Most liquids become denser when cooled and more viscous when compressed. In con- 
trast, water becomes less dense when cooled (density, or thermodynamic, anomaly) and its 
diffusivity increases upon pressurization (diffusion, or dynamic, anomaly). These anomalies, 
which disappear at high enough temperature and pressure, are not unique to water. Other 
liquids with local tetrahedral order (e.g., silica and silicon) also exhibit thermodynamic and 
dynamic anomalies A possible explanation of these anomalies is the tendency of these 
substances to form local open structures not present in simple liquids. However, establish- 
ing a precise and quantitative link between the microscopic structure and the dynamic and 
thermodynamic anomalies of tetrahedral liquids has proved elusive until recently. 

Errington and Debenedetti P| (ED) studied the relation between microscopic structure 
and the anornalies of liquid water by introducing two simple metrics: a translational order 
parameter t Q , quantifying the tendency of molecule pairs to adopt preferential separations, 
and an orientational order parameter q , quantifying the extent to which a molecule 
and its four nearest neighbors adopt a tetrahedral local structure, as in the case of hexagonal 
ice. A useful way of investigating structural order in liquids is to map state points into the 
t — q plane. Such a representation was introduced by Torquato and coworkers , who first 
applied it to sphere packings and referred to it as an order map. ED used the order map 
to investigate structural order in water j^f. Because of the distinctive features discovered in 
that study, in what follows we refer to water-like order maps as the ED order map. Using 
molecular dynamics simulation of the SPC/E [6j model, ED found that the state points 
accessible to liquid water define a two-dimensional region in the t — q plane, meaning that in 
general t and q are independently variable in liquid water (i.e., equilibrium state paths exist 
along which one order metric varies while the other does not). They also found a dome- 
shaped region in the (T, p) plane within which isothermal compression leads to a decrease in 
t and q. This decrease in order upon compression constitutes a structural anomaly: simple 
liquids, in contrast, always become more ordered upon compression. ED further found 
that dynamic and thermodynamic anomalies define nested domes in the (T, p) plane: the 
structural anomalies dome contains the dynamic anomalies dome, which in turn contains 
the thermodynamic anomalies dome. This means that whenever the thermal expansion 
coefficient is negative, the diffusivity must necessarily increase upon isothermal compression. 
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ED showed that all state points exhibiting structural, dynamic or thermodynamic anomalies 
define a line on the {t, q) plane, meaning that when water exhibits anomalous behavior, its 
translational and orientational order metrics become strictly coupled. This is clear evidence 
of the relationship between structure and water anomalies. Shell et al. subsequently found 
qualitatively similar behavior in molten silica's order map [7J. However, in the case of silica, 
it was found that state points corresponding to anomalous behavior define a narrow stripe in 
the (t, q) plane instead of a strict line. Furthermore, unlike in water, the region of dynamic 
anomalies was found to contain that of structural anomalies. 

For simple spherically-symmetric liquids, including hard spheres 0] and Lennard- 
Jones jsl the order map was found to be a positively-sloped line in the (t, q) plane, indicating 
that translational and orientational order are always strictly and positively correlated. In 
this case, of course, the appropriate metric for orientational order does not measure tetra- 
hedrality; rather, the bond-orientational order parameter introduced by Steinhardt et al. 
was used. An important result from these studies is the fact that the order map for the 
Lennard- Jones system above its critical density is identical to that of hard spheres. Further- 
more, in these simple systems that do not exhibit thermodynamic or dynamic anomalies, 
compression always leads to an increase in the order metrics. 

In 1970 Hemmer and Stell 10| showed that in fiuids interacting via pairwise-additive. 




spherically-symmetric potentials consisting of a hard core plus an attractive tail, softening 
of the repulsive core can produce additional phase transitions. This pioneering studv elicited 
^co^ide.aUe body of wo* o„ =o-ca.,ed oo.e-.ftened po.en.ia. yl^Qfly 
This generic term denotes continuous potentials with infiections in the repulsive 
discontinuous potentials with the core softened by shoulders or ramps jl^ ll^ 
17 1, or lattice models with nearest-neighbor attraction and next-nearest neighbor 
repulsion 18[ . It is now well-established that such potentials can generate water-like density 



and diffusion anomalies 



veii-estaDnsnea tnat sucn potentials 

00,0,0, 0,0,00,0. 



This important finding implies 



that strong orientational interactions, such as those that exist in water and silica, are not a 
necessary condition for a liquid to have thermodynamic and dynamic anomalies. 

The above discussion implies the existence of two well-defined classes of liquids: simple 
and water-like. The former which interact via spherically-symmetric non-softened poten- 
tials, do not exhibit thermodynamic nor dynamic anomalies, and their order map is a line. 
In water-like liquids, interactions are orient at ion- dependent; these liquids exhibit dynamic 
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and thermodynamic anomalies, and their order map is in general two-dimensional but be- 
comes linear (or quasi-linear) when the liquid exhibits structural, dynamic or thermodynamic 
anomalies. Intermediate between these well-defined extremes is the class of core-softened 
liquids, which interact via spherically-symmetric potentials but can also exhibit water-like 
thermodynamic and dynamic anomalies. 

Two questions arise naturally from this emerging taxonomy of liquid behavior. First, is 
structural order in core-softened fiuids hard-sphere or water-like ? Second, is it possible to 
seamlessly connect the range of liquid behavior from hard spheres to water-like by a simple 
and common potential,^ simply by changing a physical parameter ? 

In a recent study J^j we addressed the first question. We showed that a core-softened po- 
tential with two characteristic length scales not only can gives rise to water-like diffusive and 
density anomalies, but also to an ED water-like order map. This implies that orientational 
interactions are not necessary in order for a liquid to have structural anomalies. In this work 
we address the second question. Specifically, we use the ratio of characteristic length scales 
as a control parameter to investigate the evolution of dynamic, thermodynamic and struc- 
tural anomalies. In this manner we show that the family of tunable spherically-symmetric 
potentials so generated evolves continuously between water-like and hard sphere behavior. 
To our knowledge this is the first time that essential aspects of the wide range of liquid 
behavior encompassed by hard spheres and tetrahedrally-coordinated network-formers can 
be systematically traversed by varying a single control parameter. 

This paper is structured as follows. Sections II, III, and IV provide details on the interac- 
tion potential, simulation method, and order parameters, respectively. Results are presented 
in Section V. Conclusions and some suggestions for future work are provided in Section VI. 



II. RAMP POTENTIAL 



We perform discrete MD simulations to study the equation of state, diffusion coefficient 
and structural order as measured by the ED order map, for a fluid whose particles interact via 
a pairwise-additive, spherically-symmetric potential that gives rise to both thermodynamic 
and dynamic water-like anomalies. The model was introduced by Jagla the potential 
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energy U{r) between a pair of particles separated by a distance r is given by (see Fig. QJ. 



oo r < ao 

f/(r) = <( Uiiai - r)/ai ao < r < (1) 
r > (Ji 

The shorter distance ctq corresponds to the hard core, and the longer distance ai characterizes 
a softer repulsion range that can be overcome at high pressure. Because of its shape, it is 
called the ramp potential. The constant slope of the ramp potential for < r < ai 
keeps the force between particles / constant, so the product of separation and force rf 
will decrease when the separation r decreases. This satisfies the mathematical meaning of 
core-softening and under these conditions the thermodynamic (density) anomaly can be 
qualitatively explained by invoking the virial theorem jisl ]. 

Of interest is the ratio between the two characteristic length scales, ctq and ai, 

A = ao/ai, (2) 

which can vary between and 1. Ref. investigated the one-scale (A = 0) and two- 
scale (A = 4/7) ramp potentials. Here we investigate the full range < A < 1, with 
A = 0,2/7,4/7,5/7,6/7. 



III. MD SIMULATION 



We use discrete MD simulation; details are given in Ref. [17|. We use the NVT ensemble 
for a system composed of 850 particles with periodic boundary conditions, and we control 
the temperature with the Berendsen thermostat [2^. However, we note that we use different 
units than in Ref. the distance r, number density p, pressure P and temperature T 
are all normalized with respect to the soft core distance cxi and the potential Ui at r = 0. 
We also investigate systems with different number of particles (A^ = 1728) and confirm that 
the results do not depend on the number of particles. The simulation ranges of temperature 
and density fully cover the region where density, diffusion and structural anomalies occur. 
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IV. TRANSLATIONAL AND ORIENTATIONAL ORDER PARAMETERS 



A. Translational order parameter 

The translational order parameter 

ESQ is defined as, 

t= [ ^ \g{s)-l\ds. (3) 



Here s=rp^^^ is the radial distance r scaled by the mean interparticle distance p~^^^, p is the 
number density, g{s) is the pair correlation function, and Sc a numerical cutoff. We choose Sc 
so that it corresponds to one-half the simulation box size, and we verify that our system size 
is always large enough so that g{s) = 1 at half the box size. For a completely uncorrelated 
system, g[s) = 1, and thus t = 0. For systems with long-range order, the modulations in 
g{s) persist over large distances, causing t to grow. Between these limits, t will change as a 
consequence of the dependence of g{r) upon T and p. 



B. Orientational order parameter 



An orientational order parameter based on spherical harmonic function was introduced by 
Steinhardt et al. |pj and used in Refs. , |2l| . In this definition, all vectors connecting 

nearest neighbors (i.e., particle pairs whose separation is less than the first minimum of 
the radial distribution function) are considered. Each of these vectors, also called 'bond', 
defines an azimuthal and polar angle, and the correspondingspherical harmonic function is 
evaluated. The orientational order parameter used in Refs. jsl 0, |2l| involves the average 
of each spherical harmonic function over all bonds. 

The orientational order parameter used for water and silica in Refs. involves, first, 

the evaluation of the local tetrahedral order for each particle with respect to its four nearest 
neighbors, and then, the average of this quantity over all the molecules of the system. In 



2l| there is no such 



the definition of the orientational order parameter used in Refs. j3|,|5|,|8|, 
concept of 'local order' for an individual particle. Moreover, the number of bonds associated 
to each particle is not fixed, but instead it changes with temperature and pressure. These two 
differences led, in Ref. jl9|, to the introduction of a slightly modified version of the original 



orientational order parameter introduced by Steinhardt et al. in ref. [9]. The resu 
metric is based on the idea of a 'local order' for each particle, analogous to Refs. 



ting order 



The 



ED maps obtained with the original (global) and modified (local) definitions of orientational 



In this work we use the same order parameter introduced in Ref. We define twelve 
bonds connecting each particle with its twelve nearest neighbors. Each bond is characterized 
by its azimuthal and polar angles (6', (f) and the corresponding spherical harmonic Y£m{(^, f) 
is computed. The orientational order parameter associated with each particle i is 



Here, l^£mj^, <p) denotes the average of Y£m{0, ip) over the 12 bonds associated with particle i. 
For ^ = 6 ^3|, Qii has a large value for most crystals such as fee, hep and bcc 9]. The values 
of Q^i for each molecule in the system follow a Gaussian distribution. Q^, the averaged 
value of Q^i over all particles z, is used to characterize the local order of the system. This 
definition of order parameter is analogous to that used in water. For water, the solid at 
low pressure is hexagonal ice where each molecule has four neighbors. The orientational 
order parameter is maximum in the ice configuration and decreases as the system becomes 
less ice-like. For the ramp potential, the solid phase at low pressure has a fee structure 
where each particle has twelve nearest neighbors. Qq has a maximum value in the fee lattice 
{Q^i'^ = 0.574) and decreases as the system becomes less correlated (for uncorrelated systems, 
Qe = l/v^ = 0.289). 

V. RESULTS AND DISCUSSION 

A. Structural anomalies and order map 

Since the pair correlation function g{r) is used to compute the translational order pa- 
rameter t (Eq. 121), we first discuss the effect of density on g{r). Fig. |21 shows the effects of 
compression on g{r) at low temperature, T = 0.04, for the various values of A considered 
in this study. In all cases, there is no inner peak at r = ao for p = 1 and 1.21, and only 
the outer peak at r = ai is present at these densities. The inner peak at r = do, which is 
broad and of modest height at p = 1.66, becomes sharper and more pronounced upon fur- 
ther compression. Interestingly, structural changes brought about by compression become 
progressively longer-ranged as A increases. Thus, for A = and 2/7, the major changes in 



order are qualitatively similar. 
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g{r) involve the development of structure at length scales < o"i associated with the growth 
of the inner peak at r = a^. However, for A = 4/7,5/7 and 6/7, structural changes upon 
compression occur at distances larger than ai. In particular, for A = 6/7, the effects of 
compression are clearly discernible at r = 3cri. 

These effects of density on the pair correlation function underlie the evolution of t upon 
compression, shown in Figs. infal)-(el). Consider for example the T = 0.04 isotherm when 
A = 2/7. It can be seen that t displays a non-monotonic dependence on density: it increases 
upon compression at low densities, 1.0 < p < 1.22, decreases over the intermediate density 
range 1.22 < p < 1.76, and increases again at high densities, p > 1.76. The initial increase at 
low densities is associated with the growth of g{cTi)- The emergence of structure associated 
with the inner (hard) core causes t to decrease at intermediate densities because the initial, 
modest growth of (7 at r ~ o"o causes \g — l\ to decrease with respect to its low-density value of 
1 (see Eq. El). Upon further compression, the growth of g{<7o) above 1 eventually contributes 
additional area to the integral of |(7 — 1|, causing t to increase. This qualitative behavior of 
t is similar for A = 0,2/7 and 4/7, and is more pronounced at low T. For A = 6/7, close 
to the hard sphere limit, the pronounced growth of the inner peak upon compression gives 
rise to a monotonic density dependence of t. The case A = 5/7 is clearly transitional, with 
non-monotonic behavior at low temperature changing to hard-sphere-like monotonic growth 
of t upon compression at high temperature. 

Orientational order, as measured by Qq, shows a pronounced dependence on A, illustrated 
in Figs.|nfa2)-(e2). When A = (no hard core), Qq increases monotonically with density for 
all T. When X = 2/7, Qq begins to exhibit non-monotonic behavior upon compression. For 
this particular value of A the trend is very mild, and is best described as a virtual insensitivity 
of Qq to compression, except for an initial increase at low enough densities. For A = 4/7 and 
5/7, orientational order exhibits a marked non-monotonic dependence on density, especially 
at low temperatures. When coupled with the corresponding behavior of t, this corresponds 
to a water-like structural anomaly, whereby both order metrics decrease upon isothermal 
compression. When A = 6/7, which is close to the hard sphere value (A = 1), orientational 
order increases monotonically upon compression. Thus, there exists a narrow interval of A 
within which the ramp fluid shows water-like structural order, whereas in the pure ramp 
(A = 0) and quasi-hard-sphere limits (A ~ 1) Qq behaves conventionally upon compression. 

Cross-plotting the order metrics against each other generates the order map, whose evo- 
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lution as a function of A is depicted in Fig. HI For all values of A except 6/7, state points 
fall on a two-dimensional region, signifying that t and Qq can be varied independently. As 
is the case for silica and water 0,13]; we find, for all values of A, an inaccessible region 
where no liquid state points can be found. In the pure ramp (A = 0) case, the pronounced 
non-monotonic dependence of t on density gives rise to isotherms with well-characterized 
t-minima, the locus of which defines the boundary between the accessible and inaccessible 
regions of the order map. For A = 2/7, the barely discernible non-monotonic dependence 
of Qq on density gives rise to loops along isotherms. The non-monotonic behavior of Qq is 
fully developed for A = 4/7. This gives rise to an order map with states corresponding to 
structural anomalies lying on a narrow stripe of the order map adjacent to the boundary 
between the accessible and inaccessible regions. This behavior is strikingly analogous to 
that of water. The insensitivity of structural order to temperature, a distinguishing feature 
of hard spheres, can be clearly seen in Fig. HJ^e) by the virtual collapse of all isotherms in 
the A = 6/7 case. The transition from water-like to hard sphere order map occurs in the 
narrow interval 4/7 < A < 6/7. In particular, for A = 5/7, there is a clear evolution from 
water-like low-T behavior (T < 0.07) to hard-sphere-like high-T behavior (T > 0.07). 



B. Thermodynamic, dynamic, and structural anomalies 

We now discuss the regions of thephase diagram where structural, dynamic and thermo- 
dynamic anomalies occur. In water [2], structural, dynamic, and thermodynamic anomalies 
occur as nested domes in the (T, p) or (T, P) planes. Structural anomalies define the outer 
dome, within which isothermal compression results in a decrease of both translational and 
orientational order. Dynamic anomalies define an intermediate dome, lying entirely within 
the structural anomalies dome, and within which isothermal compression leads to an in- 
crease in the diffusion coefficient. Thermodynamic anomalies define the innermost dome, 
within which water expands when cooled isobarically. In silica [3] , dynamic anomalies define 
the outer dome, structural anomalies the intermediate dome, and thermodynamic anomalies 
define the inner dome. Thus, in both cases negative thermal expansion implies also diffusive 
and structural anomalies, but in silica diffusive anomalies occur over a broader range of 
densities and temperatures than structural anomalies, the opposite being true in water. 

Fig. shows the loci of dynamic and thermodynamic anomalies for three values of A. 
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Similar to water and silica, in ramp fluids thermodynamic anomalies occur over a narrower 
temperature and density range than dynamic anomalies. In other words, if a ramp fluid is at 
a state point where it expands when cooled isobarically, its diffusion coefficient necessarily 
increases upon isothermal compression. It can be seen that upon increasing A, the range 
of temperatures where anomalies occurs shrinks, whereas the upper limit of density (or 
pressure) where anomalies can occur increases. The shrinking of the temperature range 
where anomalies occur follows from the fact that increasing A makes the ffuid progressively 
hard-sphere-like, and there are no anomalies in a hard sphere ffuid. 

Fig. ini shows the relationship between the loci of dynamic, thermodynamic and structural 
anomalies. In water, the low-density and high-density branches of the dome of structural 
anomalies correspond to tetrahedrality maxima and translational order minima, respectively. 
For the pure ramp case (A = 0), the orientational order increases monotonically with den- 
sity over the range of temperatures explored here. Accordingly, as seen in Fig. IH^a), the 
dynamic and thermodynamic anomalies domes are bounded by loci of translational order 
extrema (maxima: line C; minima: line A). Between lines C and A, compression leads to a 
decrease in translational order. For A = 4/7 and 5/7, the locus of orientational order maxima 
(B) provides a low-density bound to the existence of thermodynamic and dynamic anoma- 
lies. Thus, for these two values of A, ramp ffuids exhibit a water-like cascade of anomalies 
(structural, dynamic, thermodynamic). For A = 2/7, maxima are barely discernible and 
the two order metrics are only weakly coupled to each other. Accordingly, the locus of weak 
orientational order maxima (B) is not a relevant indicator of dynamic or thermodynamic 
anomalies. 



VI. CONCLUSION 



In this work we have investigated thermodynamic, dynamic and structural anomalies in 
ramp potential ffuids, as a function of the ratio A of length scales corresponding to the inner 
hard core, and to the outer edge of the ramp. We ffnd that thermodynamic and dynamic 
anomalies exist for A = 0, 2/7, 4/7 and 5/7, but not for = 6/7. As in water and silica, the 
loci of anomalies form nested domes in the (T, p) plane, inside which the thermal expansion 
coefficient is negative (inner dome) and the diffusivity increases upon compression (outer 
dome). The limit A = 1 corresponds to hard spheres, and the absence of anomalies for A= 
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6/7 indicates approach to hard sphere behavior. The order map of this family of ramp fluids 
is water-like at A = 4/7 and hard-sphere-like at A = 6/7. Thus, by varying the ratio of 
characteristic length scales, the family of ramp potentials spans the range of liquid behavior 
from hard spheres to water-like. 

These findings show that orientational interactions are not necessary for the existence of 
thermodynamic, dynamic, or structural anomalies. Instead, water-like behavior apparently 
emerges in this spherically-symmetric family of fluids through the existence of two competing 
length scales. Although thermodynamic and dynamic anomalies exist almost over the entire 
range of the control parameter, the combination of thermodynamic and dynamic anomalies 
plus a water-like order map occurs over a narrow range of A. It is interesting to note that a 
distinguishing feature of water is the fact that the ratio of radial distances to the first and 
second peaks of the oxygen-oxygen pair correlation function is not ~ 1/2, as in simple liquids, 
but ~ 0.6. This is close to 0.571 (A = 4/7, the ratio of to (Xi that gives rise to water- 
like structural, dynamic and thermodynamic anomalies). In water, isothermal compression 
pushes molecules from the second shell towards the first shell, gradually filling the interstitial 
space j^l- Likewise, in the ramp potential, isothermal compression pushes molecules from 
the soft core (o"i) to the hard core (cq). Further work is needed to establish whether a ratio 
of competing length scales close to 0.6 is generally associated with water-like anomalies in 
other core-softened potentials, for example linear combinations of Gaussian [23] potentials. 
In this work we used the terminology water-like to denote structural, diffusion, and density 
anomalies. The increase in water's isothermal compressibility upon isobaric cooling, another 
of this liquid's canonical anomalies, is also trivially captured by the ramp potential, because 
thermodynamic consistency arguments j2^ mandate that the compressibility increase upon 
cooling whenever there exists a negatively-sloped locus of density maxima in the (P, T) 
plane. 

The .a.p pote„t,al, wKe. supplemented by explicit HQ o. .ean-field att.aet.o,. Q, 

gives rise to liquid-liquid immiscibility and a critical point distinct from the one associ- 
ated with the vapor-liquid transition. A liquid-liquid transition has been observed ex- 
perimentally in phosphorus j3; 27 1, n-butanol 28] and triphenyl phosphite 2^, and 
strong experimental evidence consistent with liquid-liquid immiscibility also exists for wa- 



ter 



ter j22 



30 



31 
2S, 



32 



M 



35l |. carbon j36| and wa- 



Computer simulations of silicon j33[, silica 

Q, fl Q, Q 

also indicate the presence of a liquid-liquid transition. A 
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systematic study of the effects of A and the ratio of characteristic energies {Ui and the at- 
tractive well depth) on the existence of a liquid-liquid transition, the positive or negative 
slope of the line of first-order liquid-liquid transitions in the (P, T) plane, and the relation- 
ship, if any 16j, between the liquid-liquid transition and density anomalies, would shed 



important new light on the phenomenon of liquid polyamorphism |l . l44l |45 1 . 

It is generally accepted that strong orientation-dependent interactions underlie many 
of the distinctive properties of associating, network-forming liquids such as water. Atomic 
liquids, on the other hand, exhibit simpler behavior, and in particular do not show structural, 
thermodynamic, or dynamic anomalies of the type discussed here. In this work we have 
shown that key properties of these apparently distinct categories of liquids can be bridged 
systematically by varying the ratio of two length scales in a family of spherically-symmetric 
potentials in which orientation-dependent interactions are absent by design. What other 
spherically-symmetric potentials, in addition to those possessing competing length scales, 
may give rise to water-like anomalies, is among the interesting questions arising from this 
study that we will pursue in future work. 
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FIG. 1: The middle figure shows the ramp potential with two characteristic length scales. (Tq 
corresponds to the hard core, ai characterizes the onset of soft repulsion. When A = (left figure) 

wc have a pure ramp potential (no hard core). When A = 1 (right figure) wc have a hard sphere 
potential. Ac ~ 0.6 is the ratio near which the system exhibits water-like structural, dynamic and 
thermodynamic behavior. 
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FIG. 2: Radial distribution function g{r) at T = 0.04 for different A = ao/ai values. The 
arrows indicate the direction of increasing density. The density values are 1.0, 1.21, 1.66, 2.19 for 
A = 2/7,4/7,5/7 and 1.0, 1.21, 1.66, 1.86 for A = 0,6/7. The distance r is normalized by cri, the 
soft core length. The curves are shifted vertically for clarity. 
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FIG. 3: The upper panels (al)-(el) show the density-dependence of the translational order parame- 
ter t for different A values. The solid lines are polynomial fits to the data. The lower panels (a2)-(e2) 
show the density-dependence of the orientational order parameter Qq. In each panel, the different 
curves correspond to isotherms (top to bottom) T = 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.10. 
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FIG. 4: Order map of the ramp fluid for different A. For each A, the eight isotherms (right to left) 
correspond to T = 0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.10, and the arrow indicates the direction 
of increasing density. 
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FIG. 5: Loci of thermodynamic and dynamic anomalies for ramp fluids with different A values. 
The region of diffusion anomalies is defined by the loci of diffusion minima and maxima (DM) inside 
which the diffusivity increases upon isothermal compression. The thermodynamically anomalous 
region is defined by locus of temperatures of maximum density (TMD), inside of which the density 
increases when the system is heated at constant pressure. 
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FIG. 6: Relationship between structural order and the density and diffusion anomalies in the p — T 
plane. For A = and 2/7, the domes of dynamic and thermodynamic anomalies are bounded by 
loci of t maxima (C) and minima (A), between which isothermal compression cause a decrease in 
translational order. For A = 4/7 and 5/7, the domes of dynamic and thermodynamic anomalies 
are bounded by loci of Qq maxima (B) and t minima (A), between which isothermal compression 
cause a decrease in both translational and orientational order (structural anomaly) . This cascade 
of anomalies is characteristic of water. 
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